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Abstract 

We have evaluated in closed analytical form the contribution of the three-loop 
non-planar 'triple-cross' diagrams contributing to the electron (g-2) in QED; its value, 
omitting the already known infrared divergent part, is 

a e (3 - cross) = ^(3) - — C(5) - — tt 4 + — ( a 4 + —hi 4 2 
K 1 2 sv ; 12 sv ; 135 3 V 24 

14 2, 2„ l^ox 23 2i « 47 2 113 

H 7r 2 ln 2 2 C( 3 ) H 7r In 2 tt 2 . 

9 3 w 3 9 48 

This completes the analytical evaluation of the (g-2) at order a 3 , giving 

ae(3 - lo op) = |A(3) - f «B) + % [(-. + ^ 2) - ^ 2 

239 4 139 A/oN 298 2l „ 17101 , 282591 

tt 4 H C(3) 7r In 2 H tt 2 H } 

2160 18 w 9 810 5184 J 

= ^-J (1.181241456...) . 
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Following the work of Ref. [1] we have completed the evaluation in close analytical 
form of the contribution to the electron anomaly in three-loop QED due to the triple- 
cross graphs depicted in Fig. 1). 
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Fig.l. The "triple-cross" graphs. 



The results are 



a e (3 — cross; la) = 



f C(5) - ^C(3) 



53 
2160 
4165 



7T 4 +4 



a 4 + -In 4 2 



24 



_l^ C ( 3 )_I 7r 2 ln2+ 

144 sw 6 2592 

1.285068495... + ^lnA , 



n 2 - 



515 1 



864 



In A 



(1) 



a e (3 - cross; 16) = - — C(5) + -tt 2 C(3) + — tt 4 - -n 2 ln 2 2 + — C(3) + -vr 2 ln2 
v ' ' 24 sv ' 36 SV ' 1080 6 144 sv ' 9 

1951 2 493 
~ ~648 n ~ 864 
= - 0.878968171... , 

(2) 

a e (3 - cross; 1c) = A C(5 ) - ^(3) - ^U 4 + | (a, + ^ln 4 2 ) + f A 2 2 

97 WoN 20 2l „ 1043 2 1 

H C(3) H 7r 2 ln2 n 2 

12 SV ; 9 432 48 

= - 0.026799490... , 

(3) 
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from which 



a e (3 — cross) = 2a e (3 — cross; la) + 2a e (3 — cross; 16) + a e (3 — cross; lc) 

(4) 



14 2l 2. 23 2l „ 47 2 113 , , 

= 0.785401156... + In A . 
As it is customary in this kind of calculations, A is the regularizing photon mass (in 

oo 

electron mass units), ((p) = £ l/n p , C(2) = tt 2 /6, a 4 = — -j . 

n=i n=1 ^ n 

The above results are in excellent agreement with the numerical results 



a e (3 - cross; la; Ref.2) = 1.291(7) , 
a e (3 - cross; lb; Ref.2) = -0.882(10) , 
a e (3 - cross; lc; Ref.2) = -0.021(100) , 
a e (3 - cross; Ref.S) = 0.785419(40) , 

where the In A's are understood. As the previous graphs were the last graphs for which 
the analytical value of the anomaly was still missing, on account the previously known 
results [4,5,6,7] the complete analytical expression of the anomaly in three-loop QED 
can now be written as 



/« i \ /<*\ 3 f 83 2>/0 , 215, 100 |Y l l4 „\ 
,(3-lopp)=(-) ( ^ 2 C(3)-^C(5) + — [(a 4 + -ln 4 2j 



a4 + i_l n 4 2 ^ -^7r 2 ln 2 2 



239 4 139 298 2l rt 17101 2 28259] 

7T 4 H C 3 vr 2 hi2 H yr 2 H ^ 

2160 18 w 9 810 5184 J 



= (-) 3 (L181241456...) . 

(5) 

By using the best numerical value of a e (4 — loop) = —1.557(70) , Ref.[8], and 

1/a = 137.0359979(32) , 

Ref.[9], one finds 

a e (th) = 1159652201.2(2. 1)(27.1) x 10" 12 , (6) 
to be compared with the experimental value, Ref.[10], 

a e (exp) = 1159652188.4(4.3) x 10" 12 ; 
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conversely, by using Ref.flO] as an input, one obtains 

l/a(a e ) = 137.03599941(56) 



(7) 



If p is the electron momentum and A the momentum transfer of some vertex am- 
plitude, the corresponding (g-2) contribution is extracted by keeping p on the mass 
shell, expanding the vertex up to first order in A, multiplying it by a suitable spinor 
projection operator and then performing the appropriate traces [11]. The (g-2) is then 
expressed as the sum a few hundred terms all of the kind 

/N 
d 4 k l( i%d%- , (8) 

each term having its own specific numerator N and denominator D. The numerators 
N are in general simple monomials in the scalar products of p and the loop momenta 
ki, such as for instance (p ■ ki), (ki ■ kj) etc, up to (p ■ k 2 ) 4 , while the denominators D 
have in general the form 

D = D^D^D^D^D^D^D^D^ , 

the exponents rii ranging from to 2 (0 means that the denominator is absent; the 
powers 2 arise from the expansion in A; in each term there are at most two denominators 
with rii =2); in m e = 1 units the denominators, see Fig. 2), are defined as 

Di =(p - h) 2 + 1 - ie , D 2 = (p-k 1 - k 2 ) 2 + 1 - ie , 

#3 =(p -h-k 2 - k 3 ) 2 + l-ie, D 4 = (p-k 2 - k 3 ) 2 + 1 - ie , (9) 

D$ =(p — ks) 2 + 1 — ie , Dq = k\ — ie , Dj = k 2 — ie , D$ = k 2 — ie . 

Besides the 8 denominators, there are in the problem 9 linearly independent scalar 
products p ■ ki, ki ■ kj, etc; 8 independent identities between the scalar products and the 
denominators are easily written, such as for instance 

p-fci=^(£> 6 -£>i) , 

I (10) 
fci • k 3 =- {-D 2 + D 3 -D 4 + D 7 )-p-k 2 . 

The identities can be used to express the scalar products in the numerator as a com- 
bination of the denominators, so obtaining a combination of integrands with a smaller 
number of denominators and therefore simpler from the point of view of the analytical 
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integration. It is to observed however that, as there are 9 scalar products and only 8 
denominators, one scalar product (which we choose to be p- k 2 ) must remain anyhow in 
the numerator; furthermore, the simpler terms obtained in that way can be in general 
non convergent when taken separately . 
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Fig. 2. The notation for the 3-loop integrations. 



In Ref. [1] we gave the analytical value of the simplest term 



J = ^4r- / d 4 hd 4 k 2 d 4 k 3 



7T' 



— - = 4tt 2 In 2 2 - -7T 4 . (11) 



With a straightforward extension of the results of Ref. [1] one can obtain, for instance 

p-k 2 _ 1 



Jj = / d 4 k 1 d 4 k 2 d 4 k 3 



7T 



n^n^o^ ^)"K«3). (12) 



The extension to all the other terms appearing in the expression giving the (g-2) is in 
principle straightforward but in fact rather wearisome (it implies, among other things, to 
build, for each tensor in the loop momenta, the scalar amplitudes free from kinematical 
singularities for which unsubtracted dispersion relations can be written). It was found 
convenient to work out another approach, based on the integration by part identity 
method first developed in Ref. [12]. In our case the method consists in writing 




X\d n ki 



d 



dk 



3,fJ- 



VL 



N 

15 



o, 



(13) 



where v\ is either the external momentum p or one of the three loop momenta ki, 
N/D stands for the ratio of all the pairs of numerators and denominators of the kind 
described above; note however that while in the original integrals the loop momenta have 
dimension 4, the identities are written for momenta having continuous n dimension, in 
which case their validity is obvious. The required n = 4 values are then recovered by 
means of a suitable limiting procedure. Needless to say, if the (g-2) from some graph 
is obtained as the n — > 4 limit of an n-dependent expression, this expression must be 
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evaluated for consistency in n dimensions from the very beginning; in particular, the 
traces of the 7-matrix must be evaluated in n-dimensions too. 

By explicitly performing the derivatives, using the identities (10) and rearranging 
the obtained terms in an identity like Eq.(13), one obtains a relation between various 
integrals which amounts to express one of them as a linear combination of the others. 
The number of all the possible identities is large (a few thousands, as there are twelve 
identities for each of the several hundreds N/D possibilities); they form a huge linear 
system, whose coefficients are in turn polynomials in n, from which the overwhelming 
majority of the occurring integrals can be expressed as a linear combination of a few 
basic integrals. As we are interested in the n — > 4 limit, we put n = 4 — 2uj and then 
deal with the limit by expanding everything in powers of u for ui — > . Fortunately, the 
expansion and the extremely lengthy processing of the various identities can be carried 
out in an almost mechanical way. 

A thorough investigation shows that all the three loop integrals (and therefore also 
the complete (g-2) expression in which we are interested) can be expressed in terms of 
the following 18 basic integrals: 



h 



TV 



n-2 



d n k ± d n k 2 d n k 3 



P ■ h 



h = 



h = 



h = 



tv 



n-2 



7T n ~ 2 ) / 



TV 



n-2 



TV 



n-2 



d n k 1 d n k 2 d n k 3 



d n k 1 d n k 2 d n k 3 



d n k 1 d n k 2 d n k 3 



d n k 1 d n k 2 d n k 3 



D±D 2 D 3 D 4 D 5 D 6 D 7 D 8 
1 



^-2) / 



d n k 1 d n k 2 d n k 3 



d n k 1 d n k 2 d n k 3 



TV 



n-2 



n-2 



TV 



d n k 1 d n k 2 d U k 3 



d H k 1 d n k 2 d n k 3 



D X D 2 D 3 D±D 7 D % 
1 

D X D 2 D^D % D % ' 
1 

D 2 D 3 D A D Q D 7 D 8 ' 
1 

D 1 D 3 D 4 D 5 D 7 D 8 ' 
1 

D x D z D b D 6 D 7 D % ' 
1 

D 2 D 4 D 5 D 6 D 7 D 8 ' 
1 



DiD 2 D 3 D 4 D 5 
1 

D 2 D 3 D 5 D 6 D 7 ' 
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'10 



—I 



TV 



n-2 



d n k 1 d n k 2 d n k 3 



hl = 



h2 = 



'13 



/l4 



!15 



In 



TV 



n-2 



TV 



n-2 



TV 



n-2 



— I 



TV 



n-2 



TV 



n-2 



TV 



n-2 



—I 



TV 



n-2 



d n k 1 d n k 2 d n k 3 
d n k x d n k 2 d n k 3 
d n hd n k 2 d n k 3 
d n k 1 d n k 2 d n k 3 

d n k 1 d n k 2 d n k 3 
d n hd n k 2 d n k 3 
d n hd n k 2 d n k 3 



D 2 D 4 D 6 D 7 D 8 ' 
V • k 2 



—i 



TV 



n-2 



d n k 1 d n k 2 d U k 3 



D^DsDj ' 
1 

D X D 3 D^D 7 ' 
1 

D X D 2 D A D 5 ' 
1 

D 3 D 5 D 6 D 7 ' 
1 

D 2 D 3 D A D 5 ' 
1 

D 3 D 4 D 7 D 8 ' 
1 

D 3 D 6 D 7 D 8 ' 
1 



Note that in the above list there is only one integral with all the 8 denominators (ii, 
the limit of which is Eq.(12) above), while there are no integrals with 7 denominators. 
With the exception of I\, all the basic integral are divergent in the n — > 4 limit; having 
put n = 4 — 2oj, the divergences show up as 1/cj singularities, so that evaluating them 
amounts to evaluate the coefficients of their expansions inw. As an example, the value 
of /is, which perhaps the easiest to obtain as it factorizes in the product of three 1-loop 
integrals, is 

/is = CM ( - At - 4t - - - 10 - 15w - 21w 2 - 28u 3 + 0(w 4 ) ) , (14) 
where C(u), defined as 

C(uj) = (7r"r(l + u;)) 3 , 
is an overall normalization factor, whose limiting value at u = is 1. 
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As an example of the kind of formulae obtained by expressing a non-trivial integral 
in terms of the basic integrals, let us give the relation 



TV 



j d 4 k 1 d 4 k 2 d 4 k 3 
+ I 2 (-3uj)+I 3 



(p-k 2 ) 2 

D 1 D 2 D z D A D b D 6 D 7 D s ™o 



13 o 

- — -co — 2o> — 6a;" 

4 4 



-h 



+ h + h (-w - 2a; 2 - 4a; 3 ) + / 6 Q - la; - 3o> 2 - 16or 

+ h (3cj + 14o> 2 + 64a; 3 ) + / 8 + J 9 (-gjj - ^ + ^ + 4uj2 + 176w 3 ^) 
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_1 _ 5 
~2 ~ 2 U ~ 



27 2 137 

OT 

2 2 



+ Jll hr 



12 



+ 



13 



14 



1 

6a; 
3 

64a; 2 
1 

~12a7 



-— + 

4a; 



139 1249 



36 



108 



-U! 



19225 
324 



-a; 



238333 
972 



-a; 



13 
128a; 
3 1 
8 ~ 6" 



383 „ 
+ — + 2a; 



128 



353 
76 



3041 



-a; 



-a; 



29 2 37 < 

— or H ar 

24 4 



15 



13 87 

T + ¥" + 

11 19 



439 2 623 
—a; 2 + ■ 



55 



73 



12u 



16 



-7^-7-7^+V + 196^ 3 
2a; 4 4 4 



17 



2^ + 4 + 



29 



50 



— a; H a;' 

6 3 



18 



64a; 2 



2 

3a; 



55 
36 



1391 

864 



-u — 



94423 2 3535507 
-a; 2 + 



5184 



7776 



LO' 



(15) 

Similar expressions hold for all the other integrals contributing to the anomaly. In all 
those expressions, the coefficient of 1\ is always finite (i.e. not singular for uj — > 0), so 
that only its u> = limit, given by J\ of Eq.(12), is required. 

Besides ii, we have already given the explicit analytical value of i"i 8 , Eq.(14). In 
order to complete the calculation of all the basic integrals, we have evaluated directly 
only a subset of them (namely J 8 , J 9 , I\q and from J 13 to I\% ). For the others, we found 
it more convenient to write a large number of identities of the form of Eq.(15), in which 
however the analytical value of the integral on the l.h.s. is already known from previous 
work; the relation is then expressing such a known result in terms of the still unknown 
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basic integrals. In the case of J , Eq.(ll), one finds for instance 



Jo = lim 



15 



-co 



'l 3 x 

— -UJ 



+ h \2 2 . 

,11 4 ,„ 



Mi 



( i 



+ i4 



16 64 9 «\ /ll 

39 o 201 -A 



5 5 \ 
A A J 

- + -uj - 5uj 2 - 20uj 3 
2 2 



1_ 1_ ; ( ; N 



■CO- 



LO 



T 1 49 11 

+ In 1 ■ 

\6cv 36 108 

/ 3_ 59 

+ 12 V 128w 3 + 256c; 2 

T ( 1 17 211 

+ il3 



1819 2 34249 -\ 

uj 2 u 6 

324 972 J 

203 15 37 53 2 

"16" 



+ 



" 16 W 



256w 

17 211 423 2 1263 , 
—— — to — — w — a; 



113 



48 16 



1 

3w 24 
1 49 

32w + 32 8" ' 16 



53 3 
— — —to + 



119 2 ^ 427^ 3 X 



3 

■a/ + 



4 17 

3uj 6 



211 423 2 1263 , N 

a; uj uj 

12 A 2 



r / 1 13 7 5 9 2 67 o\ 
+ /l6 l8^-16^ + 16-r-2 Ci; "T W ) 
T ( 3 23 69 40 587 ,\ 

+/l7 ( v -4^ + 8 ^-y-y u; --6- Ci; J 

/ 3 7 29 3799 140051 1193035 2 10388627 

+ 18 V 128w 3 + 32w 2 ~Y2u~i~ ~576~ ~ 3456 U 5184 U 7776 U 

Similarly, in Ref.[7] the following 7-denominator integral was evaluated 



1 <? Mk > b^oAd^d, = §^< 3 > - f C(5) ; 
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in terms of the basic integrals one has 




(— 7 r f 1 

+i "{-l +1 2--^ 2 - i ^) +i '{-^ + l) 




, ! 1 11 1 2 i 3 \ , ( i i 




1 1,1 




(18) 



It is to be stressed here that such relations have quite a broad generality, which extends 
to cover scalar amplitudes occurring in graphs with different topology; the l.h.s. integral 
of Eq.s(17),(18) was indeed evaluated in Ref.[7], dealing with "corner-ladder" graphs. 

By exploiting our previous work, we can write a redundant number of such re- 
lations expressing unknown basic integrals in terms of already known integrals (they 
are not written here for the sake of brevity); in so doing we obtain a redundant set 
of linear equations, which are easily solved in terms of the unknown basic integrals 
(the redundancy provides with additional consistency checks). As a side remark, let 
us observe that for each basic integral one must evaluate the first few terms (up to 7) 
of its expansion in u>, see for instance lis , so that the number of unknown quantities 
increases, but on the other hand each relation, being an identity in u> provides with 
several independent equations for the unknown. 

As a result the following integral table is established 



h = C{u) 5C(5)--tt 2 C(3) + OM 



h = C(u) 2 




3^ 2 + 10C(3) 

Fi 7 \ 

-tt 2 C(3) - — tt 4 - 82C(3) - 47r 2 ln2 + 16tt 2 + 0(u 2 ) 
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h = C(u) 
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he = C(u) 
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2u 3 Au 2 
16 



7 1 {I o 25 
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By using the above table in combination with the integration by parts identities, we 
obtain the analytical value of the required integrals; from Eq.(15) we have for instance 



(-if 



d 4 k 1 d 4 k 2 d /L k 3 



(p ■ k 2 f 



DiD 2 D 3 D 4: D 5 DqD 7 D s 



5C(5)- A 7r 4 + c(3) _ 27r 2 ln2+ 3 7r2 



(18) 
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In terms of the basic integrals, the (g-2) of the graph of Fig.lc), for instance, reads 



(<7-2)(Fig.lc))= lim 



' 1 19 
+ 12 



1637 
^72 



-uj 



, ( 13 

1 



+ 



+ / 6 | 



253 
144 
7_ 

18 ~ 

55 



823 102797 2 

o> H uj 

108 2592 

587 



979525^ 3 N 



~2Auj + 72 



( 19 25 632 

/ 1 5 1585 _ 
+ 8 V8^2 ~ 9^ + "216" ~ 1296 



3888 

9133 o 340685 ,\ 
324 1944 J 

19075 o 874721 3 X 

56983 2 , 1488295^ 3 \ 



+ / 4 



v 8o> 



5 229 N 
24 + 18^ 



216 
827 

108"^ 648 



648 
36581 



u 2 + 



■uj + 



1944 

1051253 2 743606 3 N 

a; 2 uj 6 

729 



+ h( 
+ ho 
+ hi 

+ Il2 



11 

96a; 2 " 

(-— 

\ 72a; 

(-— 

V 6a; 2 
11 



565 467 11225 

+ — TT^U + 



288a; 
_ 613 
" 216 
1 



36 
187 



3888 

1629641 2 

a; 2 - 



3888 



1068433 
432 



■uj c 



324 
2143 
18u ~ ~216~ 
525 



■uj — 



+ /l3 



11 525 6601 

512a; 3 + 1024a; 2 ~ 2304a; 
1632171647 3 \ 
373248 U J 

95 

96a; 2 
/ 35 



162 

234293 2 1054147 3 \ 

3888 23328 J 

5767 2472821 2 20840363 ,\ 

+ — "--7776-" +^1664-" ) 

577597 4981223 46026487 2 

a; H uj 

41479 ^ 62208 



144 
6601 



+ 



27648 " 41472 



■ + 



+ 



192a; 

/ , .,, _ 1051 

+ 14 V192a; 2 1152a; 



339883 2 60055 , 

uj H uj H a; 

576 432 2592 3888 

29419 220349 

H a; + 

1728 3456 



264997 2 31472155 c 

w " -TTT^l— " 



+ /i 5 ( 



1 7 

+ 



,4a; 2 
' 11 
96a; 2 ~ 

' 17_ 

24a; 2 
11 



36a; ~'~ 


432 + 


UJ + 

27 


1075 


8179 


10925 


576a; 


+ 864 


192 



3888 
■uj + 



486 

T 11664 



1934333 2 , 1580447^ 3 N 



2919977 2 49397413 ^ 
a; 2 u 6 



1429 12223 



+ 144a; 
1315 



+ 



"512a; 3 ^ 1536a^ 
177573299 3 Y 
" 93312 U ) 



216 

3107 ( 292849 
~ 2304a; 



7776 " 23328 

15117413 2 X 
432 " 7776 

206467 



146749 
uj — 



-a; z 



6912 4608 



-uj 



181740889 2 

a; 2 

124416 



(19) 



14 



Similar expressions hold for the other graphs. 

By using the above table Eq.s(2,3) are immediately obtained. 

Graph la is ultraviolet divergent and requires renormalization. From the unrenor- 
malized amplitude we obtain 



1 215 1 
a e (la;not ren.) = - — + ^jC(5) - ^Ci 3 ) 



53 4 „ 

7T 4 + 4 



2160 



a 4 + — (In 4 2 - 7r 2 ln 2 2) 



(20) 



1229, , 7 2l 3571 2 133 
l4T C<3) - r 1112 + 2692" + 864 



To obtain Eq.(l), one must account for the charge renormalization of the inserted 4th 
order vertex amplitude, to be carried out by means of a suitable subtraction constant, 
which will be indicated here as Z. For consistency with the previous analytic (g-2) 
calculations, Z must correspond to on mass-shell renormalization. It turns out therefore 
that Z is infrared divergent (this is the only part of the triple-cross (g-2) calculation 
in which infrared divergences appear); again for consistency with previous work, the 
infrared divergence is parametrized by giving to the photon a fictitious infinitesimal 
mass A. The ultraviolet divergences of the counterterm are still parametrized by means 
of the n-dimensional regularization; the result for the renormalization counterterm of 
the inserted 4th order graph reads 

z =—z- ~ ^r 2 + 2-lnA ■ 
8u 24 

In order to renormalize the contribution of graph la, one must subtract from Eq.(20) 
(0)Z where (0) is the second order (1 loop) magnetic form factor, to be evaluated 
in n-dimensions as well; as Z has a 1/uj singularity, use must be done of the value of 
F^ ; (0) evaluated up to first order in u, 

F 2 (2) (0) = l + 2^. 
By subtracting F^\q)Z from Eq.(20), Eq.(l) is recovered. 
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